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Abstract 

We study an exclusion process with 4 segments, which was recently introduced by T Banerjee, N Sarkar 
and A Basu [J. Stat. Mech. (2015) P01024]. The segments have hopping rates 1, r(< 1), 1 and r, 
respectively. In a certain parameter region, two shocks appear, which are not static but synchronized. 
We explore dynamical properties of each shock and correlation of shocks, by means of the so-called 
second-class particle. The mean-squared displacement of shocks has three diffusive regimes, and the 
asymptotic diffusion coefhcient is different from the known formula. In some time interval, it also 
exhibits sub-diffusion, being proportional to . Furthermore we introduce a correlation function 
and a crossover time, in order to quantitatively characterize the synchronization. We numerically 
estimate the dynamical exponent for the crossover time. We also revisit the 2-segment case and the 
open boundary condition for comparison. 


1 Introduction 

The exclusion process is a paradigm for non-equilibrium behaviour [1] . In the one-dimensional totally asymmetric 
simple exclusion process (TASEP), each site i (* = 1, 2,..., L) is either occupied by one particle (t^ = I) or empty 
{ji = 0), and a particle at site i hops to j -I-1 with rate I, if the target site is empty. Investigation of shocks is one 
of central issues in the TASEP n in nnainiH]. Let us recall a known result in the TASEP with open boundaries 
(shortly, open TASEP), where a particle is injected at site i = 1 with rate a, and extracted at site i = L with 
rate /3 [9]. These rates are regarded as reservoir densities a and 1 — /3, respectively. The case a = /3 < 1/2 
is called the co-existence line. Two plateaus with densities a and 1 — a co-exist in the domains 1 < x < S{t) 
and S{t) < X < L, respectively. The position S{t) of the shock (domain wall) is dynamical, and its behaviour is 
diffusive 11, *.e. 


{{S{^ - S{0)f)^ 2D{a)t, D{a) = (1) 

as t —>■ 00 . Here ( ■ )g denotes the ensemble average. This asymptotic diffusion coefficient is also true on Z 1. 

The so-called second-class particle (r^ = 2) microscopically defines the positions of shocks [5]. It behaves as a 
hole for particles, and as a particle for holes, i.e. the hops 20 —>■ 02 and 12 —)■ 21 occur between sites i and i -I-1, 
with the same rate as for 10 —>■ 01. We conhne only one second-class particle into the system in the initial state. 
It is not extracted from the system and we do not inject another second-class particle, i.e. “semi-permeable 
boundary condition” uni El m [131 • Figure0(a ,b) shows simulation results of the mean-squared displacements 
(MSD) of the second-class particle, which agree with the formula Q. We notice that finite-time effect becomes 
strong, when a approaches 1/2. The density profile averaged over a large time interval 

Pi ■= {Ti(2. - Tj))t (2) 

looks different from snapshots in simulations. Since the shock position moves evenly in the chain, pi is given by 

p^~{l-2a)^+a, (3) 

which was shown by the exact stationary state |3]. 
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Figure 1: (a) MSD of the second-class particle vs time, and (b) diffusion coefficient vs boundary rate on the 

co-existence line of the open TASEP with L = lO'^. The markers in (a) are simulation results averaged over 10^ 
runs. We used simulation data (S'(t) — S'(0))^/(2t) in 5 x 10^ < t < to plot the markers in (b). The lines in 
(a,b) correspond to eq. Q. 


On the other hand, static shocks are realized in the TASEP, by imposing attachment and detachment of 
particles in the bulk of the chain (the Langmuir kinetics [II]). Static shocks were also found in TASEPs with 
inhomogeneous hopping rates on a ring. One of the simplest cases is the Janowsky-Lebowitz (JL) model [15) : 

Pi = 1 (1 < * < A), r < 1 {i = L), (4) 

where pi denotes the hopping rate from site (A+1 := 1). Due the inhomogeneity of the bond between sites 

L and 1, the JL model exhibits a shock profile in a certain parameter region. A mean-field theory qualitatively 
explains a phase transition between shock and flat density profiles, which is still a fascinating problem [Mill]- 
In [H], another inhomogeneous TASEP was introduced: 

Pi = 1 < i <(.), r < i < = L). (5) 

We refer to this model as 2-segment TASEP. In a certain region of the parameter space (r, p) (where p is the 
global density i.e. the number of particles /A), this model also exhibits a static shock. Recently the authors of 
m investigated TASEPs with three and four parts, generalizing ([^ . Here, we mainly study a specific 4-segment 
TASEP 


Pi = 


V 2^<i< 3£), 

{i<i<2i V i < 4£ = L). 


( 6 ) 


There can exist two shocks in the 1st and 3rd segments. The positions of the shocks cannot be fixed even in the 
macroscopic level, but they are related to each other by an equation derived by the particle number conservation 
m- The purpose of this work is performing more detailed Monte Carlo simulations (in continuous time), in 
order to deeply understand this synchronization phenomenon. 

Before investigating the 4-segment TASEP, we reconsider the case of 2 segments via the second-class particle. 
As an evidence that the second-class particle indicates the shock position, we check that its spatial distribution 
becomes gaussian, corresponding to the density profile written in the error function. We also examine properties 
of the standard deviation of the shock position. Then we turn to the 4-segment TASEP. We find that the MSD 
exhibits various behaviours, depending on the time scale that we focus. The open boundary condition that we 
have already reviewed is the reference case. We quantify the interaction between the shocks by a correlation 
function. Furthermore we introduce a crossover time distinguishing between time scales of independency and 
synchronization of shocks. We also numerically estimate the dynamical exponent of the crossover time. Finally 
we give the conclusions of this work and some remarks, including possible future studies. Overall in this work, 
we use the following definition for macroscopic density profiles with mesoscopic length of the lattice 2d -|- 1 
(1 <C d <C £), which is in general different from the microscopic density profile ([^: 

, i-l-5 

p{x) = p{i/e) = ^ Ti,i2-Ti>). (7) 

i'—i—S 
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Figure 2: (a) Phase diagram and (b) fundamental diagram of the 2- and 4-segment TASEPs [HI [19]. The phase 
boundaries in (a) are given hy p = Pc, I — Pc *-e- the parabola r = 1 — 16(/9— 1/2)^. For (b), we have set £ = 10^ 
and r = \I2. To obtain each marker, we counted the number of flowing particles in a single simulation run, and 
took average over 10® <t< 10^. The line in (b) corresponds to the predictions (12), (16). 


2 2-segment TASEP 

Let us consider the model ([^. We begin with the assumption that the global density p is enough small, and each 
segment has a flat density profile 


p{x) = ai(0 < a; < 1), 02(1 < x < 2). 

The conservation laws of the number of particles and the stationary current J are written as 


2p = ai -I- a2, 

J = ai(l — ai) = ra2(l — 0 : 2 ), 


respectively. These are easy to solve m- 


1 -I- r — 4rp — R 4p — {1 + r) + R 

= - 7 ^ 7 ^ -\-) 0^2 = 


J = 


2(1-r) ’ " 2(1-r) 

r(l-2p)(i?-(l + r)(l-2p)) 


(1 - r )2 

with R = \/{\ + rY — 16rp(l — p). The density profile (|^ is realized as long as 02 < 1/2 with i.e. 

2-Vl^ 


p< Pc-= 


( 8 ) 


(9) 

( 10 ) 


( 11 ) 

( 12 ) 

(13) 


Since both ai and q ;2 are less than 1/2, the case p < pc is called LD-LD phase (LD=low density). 

When the total density p exceeds the critical density, the 2nd segment maintains the density 1/2 and a shock 
appears on the 1st segment: 


p{x) = 


Q !1 (0 < X < s), 1 — «! (s < X < 1), 
Q!2 = I (1 < X < 2). 


(14) 


We refer to this case as S-MC (shock-maximal current) phase. The shock position s is macroscopically static, 
i.e. localized. By solving the conservation laws of the number of particles and the stationary current 


2p =sai + (1 — s)(l — Qfi) -I- 0 : 2 , 
J{ai) =Q!i(l - ai) = ra 2 (l - 0 : 2 ) = r/i, 

the densities and the shock position in the 1st segment are determined as 



(15) 

(16) 


(17) 


The parameters (r, p) are inversely specified by (ai,s) in the S-MC phase. We have emphasized the current as 
a function of the density ai. 
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Figure 3: (a) Distributions of the second-class particle, (b) density prohles, (c) standard deviation a vs segment 
length I, and (d) cr vs hopping rate r. For (a) and (b) we have set the system length as£ = 200, and for (d) £ = 10^. 
Each plot marker was obtained by averaging data of a single simulation run over 10® < t < 10^ [exceptionally 
10® < t < 5 X 10^ for the cases £ > 10^ in (c)]. The curves in (a) and (b) correspond to eqs. (|I^ and (19), 


respectively, with (S) and a obtained from simulation data. Each line in (c) corresponds to {£/\Q'^) 
with ah^iQi obtained by simulations. 


4Tr72 


X cr|i=io4 


When the global density p exceeds 1 — pc, the shock position reaches site i = 1. The densities in both 
segments are flat, and larger than 1/2, i.e. the HD-HD phase (HD=high density). Figure (a) summarizes 
the three phases. In hg. (b), we check that the predicted currents (12) and (16) are realized by simulations. 
Because of the particle-hole symmetry, we restrict our consideration to p < 1/2. 

Let us investigate properties of the shock in the S-MC phase in more detail. There is a microscopic deviation 
of the density profile near the shock position. It obeys a gaussian distribution, i.e. the probability distribution 
P{S) := \ {ts{ts — I))t of the position S of the second-class particle is given as 


P{S) = 


\pli 


: exp 


2(t2 




(18) 


with (S')t ~ £s {£ ^ oo). We see good agreement of simulations to this distribution in fig.[^(a). This corresponds 
to the fact that the (microscopic) density profile ([^ in the 1 st segment is well described in terms of the error 
function erf[-] [19], see hg. |^(b): 


Pi = —s— erf 


»-((^)T+i)^ 


(19) 


We chose the values of parameters, such that corresponding densities and shock positions are given as 


(r, p) = (0.96,0.5), (0.84,0.48), (0.64,0.44), (0.36,0.38) 

(ai, s) = (0.4,0.5), (0.3,0.6), (0.2,0.7), (0.1,0.8), (20) 

respectively, according to eq. (In the 2nd segment, we expect that the density proHle, which is deviated 

from 1/2, can be well written by the exact hnite-size effect in the maximal current phase of the open TASEP 

El-) 

In hg. (c), we observe cr oc [19]. Furthermore a/VI vs r is shown in hg. (d) with different global 
densities p. So far we have not found an explicit formula of cr/-\/Z, but it seems independent of p. 


3 4-segment TASEP 


Now we turn to the 4-segment TASEP ® m- By the same argument as for the 2-segment case, the density aj 
of each segment j is hat, when p < pc. (The phase transition line p = pc is identical to that of the 2-segment case, 
hg. I (a).) From the translational invariance, we have ai = 0:3 and a 2 = a 4 . Furthermore these densities have 
the same form (11) as for the 2 -segment case [TU]. After the global density exceeds pc, the 2 nd and 4th segments 
maintain the density 02 = a 4 = 1/2, and shocks appear in the 1st and 3rd segments [IH]: with ai = as, 


{ ai (0 < X < si), 1 — «! (si < X < 1), 

i (1 < X < 2), 

as (2 < X < S 3 ), 1 — as (sg < x < 3), 
I (3 < X < 4). 


( 21 ) 
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Figure 4: (a) Kymograph probing second-class particles and macroscopic density profiles [eq. Q with d = 5] 

in the 1st and 3rd segments for {r,p,£) = (0.64,0.5,200). (b) Distribution of the second-class particles, and (c) 
density profiles for £ = 200, averaged over 10® < t < 10®. The lines in (b) are the uniform distribution over the 
range (24), and in (c) the prediction ([^ with x = il£. 


The shock positions are denoted by si and S 3 = S 3 -I- 2 [53 = 0 (resp. S 3 = 1) corresponds to the boundary 
between 2nd and 3rd (resp. 3rd and 4th) segments]. The conservation of the number of particles is written as 


4 p = siai -I- (1 - si)(l - ax) + a2 + sza^ -b (1 - S3)(l - as) -b 04. 


( 22 ) 


Solving this together with the current conservation (10 1 , we find a restriction on the shock positions |19] 


Si -b S 3 = 2 s, 


(23) 


with s The form of the density ai is the same as for the 2-segment case ( [T7| ). The current J is also 

unchanged m, see fig. Sb). One of the interesting findings in [T5] is that we cannot fix the shock positions Si 
and S 3 even in the macroscopic level, but they are synchronized by the restriction ( [2^ , see fig.|^(a). The right 
bound of Sj is 1, and the left A is given by solving 2s — A = 1: 


2(1 - 2p) 

•\/l — r 


< s, < 1. 


(24) 


The shock positions are uniformly distributed in this range except for the boundaries, see hg. (b). Therefore 
the density profile by averaging configurations over a large time interval becomes 


(p(x))t = 


{ ai (0 < a; < A), ai -b (x — X)g (A < x < 1), 

1 (1 < X < 2), 

ai (2 < X < A'), ai -b (x — X')g (A' < x < 3), 
1 (3 < X < 4), 


(25) 


where A' = A -b 2 and g = 
the prediction (251 with x 
left boundaries become 


In fig. 1^ (c), we compare simulation results of the density profile P. (H with 
= il£. We chose the values of parameters, such that the corresponding densities and 


(r, p) = (0.96,0.5), (0.84,0.48), (0.64,0.44), (0.36,0.38) 

(ai. A) = (0.4,0), (0.3,0.2), (0.2,0.4), (0.1,0.6), (26) 


respectively, according to eqs. (|T7| and (24). We expect that the discrepancies between the prediction and 


simulations are due to the finite-size effect. In the half-filling case p = 1/2, the profiles of time average in the 1st 
and 3rd segments are identical to that of the well-studied open TASEP Therefore one may naively expect 
that the 1st and 3rd segments are effectively equivalent to the co-existence line of the open TASEP; e.g. for the 
3rd segment, the 2nd and 4th segments play the role of reservoir of densities ai and 1 — ai, respectively. Let us 
observe simulation results of MSD in the 4-segment TASEP, to verify whether this guess is correct or not. 

Denote by Sj{t) + {j — 1)^ the microscopic position of the second-class particle in segment j S {1, 3} at time 
t, so 1 < Sj{t) < £. Figure]^ (a) shows the simulation results of the MSDs 




E' 


(27) 
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Figure 5: (a) MSD of the shock positions, (b) initial diffusion coefficient, and (c) exponent e (31). In (a), we 

divided M{t) by 2t, so as to see its property more easily. In (a,b), the solid, dashed and dotted lines correspond 
to J [eq. ([T^], D [eq. 0] and V [simply M(10®)/(2 x 10®) from simulations], respectively. To obtain numerical 
data in (a,c), we took averages over 10® and 10® simulation runs for t < 2 and t > 2, respectively. For each 
marker in (b), we averaged M{t)/{2t) of 10“® <t< 10“^ over 10® simulation runs. 


Thanks to the translational invariance, we omit the subscript j, and practically we obtain M(t) by averaging 
simulation data of j = 1 and j = 3. In some time interval, indeed we see that the diffusion coefficient ([^ with 
a = ai 0 provides a good estimation of the MSB’s behaviour: 

M(t) « 2D{ai)t. (28) 

Note that the observed M(t)/(2t) in this intermediate diffusive regime is slightly lower than D, which is expected 
to be a finite-size effect. 

Figure (a) indicates that the initial and asymptotic behaviours are also diffusive. The current of particles 
0 gives a good fitting curve for the diffusion coefficient in the initial diffusive regime^ see fig. [5](b): 

M{t)~2J{ai)t (t^-fO), (29) 

which implies that the second-class particles obey the symmetric random walk with jump rate J(ai) in a very 
short time scale. The same finite-time effect is observed in the open TASEP [fig. (a)], where one can prove 
that both probabilities of finding configurations 20 and 12 are given by J{a) by using matrices [TUI ITTl IT^ IT^ . 
On the other hand, for the asymptotic diffusive regime, 


M{t) ~ 2V{ai)t {t oo), (30) 

we expect the inequality T>{ai) < D{ai), see fig.[^(a). Note that the conclusion on the asymptotic behaviour is 
based on our simulations of up to < = 10®. Intuitively, no further transition is expected in longer-time simulations. 
(We considered the “asymptotic” behaviour in t ^ ^I'D, so that the shocks do not reach boundaries Sj{t) = X£,£. 
When we take the limit t ^ oo with £ fixed, the MSD converges to a constant.) The notation T>{ai) emphasizes 
that the diffusion coefficient is a function of ai, but the dependence on £ and A may also be nontrivial. 

In the time interval between the intermediate and asymptotic diffusion regimes, the MSD exhibits sub¬ 
diffusion. The bold lines (slope —1/2) are drawn in the log-log graph, fig.[^(a), which agree well with simulations. 
We also estimate its exponent e (In M (f) « e In t) by 


e = 


1 

3t/2 -h 1 


E 

t/2<t'<2t 


InM(lOF) -lnM(F) 
ffilO 


(31) 


in this transient regime. The exponent by this formula takes the minimum, which is « 1/2, see fig.[^(c). 

Let us investigate the correlation of the shocks. We expect that, in a very short time, the two shocks move 
independently, since they are far from each other. On the other hand, they are synchronized in the frame of 
a larger time scale, e.g. as hg. (a). In order to quantify (in)dependency of the two shocks, we introduce a 
correlation function 


C{t) = -((5i(t) - Sm){S3{t) - Ss{0)))^. 
Due the synchronization Si(t) + S^ft) « £s, we have 

C{t) oi M(t) (t^+oo), 


(32) 

(33) 
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Figure 6: (a,b,c) Comparison between the MSD and the correlation function in different time windows, and 

(d) crossover time vs segment length. For (a,b,c), we have set the values of the parameters as = 

(1000,0.64,0.44), and averaged over 10^, lO'^ and 10® simulation runs, respectively. In (d), each marker was 
plotted by averaging over 10® simulation runs, and the lines are htting curves ci^. 


which is observed in fig. (a). On the other hand, their initial behaviours are completely different, see hg. 
(b,c). In the vicinity of t = 0, we conjecture that C{t) increases more slowly than any power function t“(a > 0). 
Let us define the time T as 


T = inf{t >0|2C'(t) >M(f)}, 


(34) 


characterizing the crossover between the time scales of independency and synchronization. Figure [^(d) shows T 
vs segment length i. By assuming the power law T c£^, we perform fitting from the simulation data that we 
have. We write the results directly in fig.[^(d), and draw corresponding straight lines as well. Under the further 
assumption that the exponent is independent from parameters, we simply average the four obtained exponents: 
z « 1.71. 

4 Discussions 


In this work, we investigated synchronization of shocks in the 4-segment TASEP, by means of the two second-class 
particles. We found that the behaviour of the MSD M{t) of the shocks is not so simple. In the initial regime 
(very short time scale), it is diffusive and the coefficient is nothing but the formula for the particle current. In 
the intermediate diffusive regime, the coefficient D known in the open TASEP provides a good estimation. This 
fact indicates that, up to the intermediate diffusion regime, the shocks’ motions are determined by the values of 
low and high densities, and independent from details of boundary conditions. Via sub-diffusive regime, the MSD 
achieves the asymptotic diffusive regime^where the diffusion coefficient is different from D. From the log-log 
graph and the numerical estimation, fig. N (a,c), it seems that M{t) oc 'Jt in the sub-diffusive regime. In some 
random walks, similar changes of the diffusivity have been found EOllll]- In our case, the regime change is 
spontaneously induced by the particle number conservation, without directly imposing any interaction between 
the two second-class particles in dehning the model. 

We also investigated the correlation function of the shocks, and the crossover time between independency and 


synchronization of the shocks. The correlation function C{t) (32) becomes identical to M(t) as t ^ -|-oo. In the 
vicinity of t = 0, C{t) increases very slowly, whereas M{t) is linear. We defined the crossover time T as the time 
when the ratio C/M exceeds 1/2. We estimated the dynamical exponent z for T by using our simulation data, 
which was found between the KPZ z = 3/2 [22] and the normal diffusion z = 2. The exponent z as well as V 
and e should be more precisely estimated in larger systems with longer simulation time and by more simulation 
runs. 

The generalization to the 2n-segment case is straightforward: the jth segment has the rate 1 (resp. r) when 
j is odd (resp. even). When p < pc, the density profile p{x) (j — 1 < x < j) of jth segment is give as p{x) = ai 
for odd J, and p{x) = a2 for even /, with the same forms as in the 2-segment case When the global density 
Pc < p < £ — Pcj the density profiles are flat with density 1/2 in the segments with even numbers. On the 
other hand, n shocks appear in segments j = 1, 3,..., 2n — 1. Denoting their positions by Sj + j — 1, we find 
Si + S 3 S 2 n-i = ns, from the conservation of the number of particles. As an analogue to the n = 2 case, we 

expect that they are not static but synchronized. Since only one equation governs the synchronization of the n 
shocks, we naturally expect that the asymptotic diffusion constant I?„(a) enjoys 'Di{a) < I? 2 (a) < ^^{a) < ■ ■ ■ 
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with = 0 and I? 2 (c«) = 'D{q). A further intuitive conjecture is that the sub-diffusive regime vanishes as 
n —>■ -1-00, and hm„_>oo T^nioi) = D{a). 

It is known that two shocks are synchronized in the model [23j. A simple generalization of the JL model 
|24| . e.g. Pi = r{i G {£, 2£}), l(i ^ {£,2£}) with L = 2£, also exhibits the same type of synchronization. One of 
important questions is whether the MSDs of shocks in these models behave like fig. (a), and if yes, whether 
the asymptotic diffusion coefficient is the same as for the 4-segment TASEP. Note that the positions of the two 
synchronized shocks are, in general, far from each other in these models as well as the 4-segment TASEP. (See 
e.g. [25l [26l [271 [28] for other types of synchronization.) While the motions of the second-class particles are locally 
defined, the shocks move as if they knew each other’s position. We believe that this viewpoint gives hints to 
study self-organization phenomena, e.g. in biological cells, as the exclusion process is one of basic models also in 
biophysics [75]. Application to traffic flows would be also an interesting problem. The local inhomogeneities in 
the JL and generalized JL models are very similar to traffic lights, and combinations of segments that we studied 
here evoke different limit speeds of cars. 
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